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ABSTRACT 

Jupiter Trojans are thought to be survivors of a much larger population of 
planetesimals that existed in the planetary region when planets formed. They 
can provide important constraints on the mass and properties of the planetesimal 
disk, and its dispersal during planet migration. Here we tested a possibility that 
the Trojans were captured during the early dynamical instability among the outer 
planets (aka the Nice model), when the semimajor axis of Jupiter was changing 
as a result of scattering encounters with an ice giant. The capture occurs in this 
model when Jupiter's orbit and its Lagrange points become radially displaced 
in a scattering event and fall into a region populated by planetesimals (that 
previously evolved from their natal transplanetary disk to ~5 AU during the 
instability). Our numerical simulations of the new capture model, hereafter jump 
capture, satisfactorily reproduce the orbital distribution of the Trojans and their 
total mass. The jump capture is potentially capable of explaining the observed 
asymmetry in the number of leading and trailing Trojans. We find that the 
capture probability is (6-8) x 10^'' for each particle in the original transplanetary 
disk, implying that the disk contained (3-4) x 10^ planetesimals with absolute 
magnitude H < 9 (corresponding to diameter D = 80 km for a 7% albedo). The 
disk mass inferred from this work, Mjisk ~ 14-28 MEarth? is consistent with the 
mass deduced from recent dynamical simulations of the planetary instability. 



1. Introduction 

Jupiter Trojans are populations of small bodies with orbits similar to that of Jupiter. 
They clump around two equilibrium points of the three-body problem, known as the La- 
grange L4 and L5 points, with semimajor axes a ~ ajup, eccentricities e < 0.15, inclinations 
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i < 35°, and AA = A — Ajup ~ ±60°, where A is the mean longitude and index "Jup" denotes 
Jupiter. The angle AA oscillates with a period of ~150 yr and full libration amplitude up 
to ~ 70° (see Marzari et al. 2002 for a review). The total population of Jupiter Trojans 
is estimated to be ^10% of the main asteroid belt (e.g., Jewitt et al. 2000, Nakamura & 
Yoshida 2008). 

Jupiter Trojans are thought to have been captured from a much larger population of 
small bodies (planetesimals) that existed in the planetary region (~ 5-30 AU) when the 
giant planets formed. Previous theories suggested that the Trojans were captured during 
the early stages of Jupiter's growth (Marzari & Scholl 1998, Fleming & Hamilton 2000), by 
collisions (Shoemaker et al. 1989), effects of nebular gas (Yoder 1979, Peale 1993, Kary & 
Lissauer 1995, Kortenkamp & Hamilton 2001), etc. These theories imply that the inclination 
distribution of the Trojans should be relatively narrow with most orbits having i < 10° 
(Marzari et al. 2002). By contrast, observations show a wide inchnation distribution of 
Jupiter Trojans with inclinations up to ~ 35°. Attempts to explain the large inclinations 
of Trojans by exciting orbits after capture have been unsuccessful, because passing secular 
resonances and other dynamical effects (e.g.. Gomes 1998, Petit et al. 1999, Marzari & Scholl 
2000) are not strong enough. 

Morbidelli et al. (2005, hereafter M05) proposed that Jupiter Trojans were trapped in 
orbits at L4 and L5 by chaotic capture. Chaotic capture occurs when Jupiter and Saturn 
pass, during their orbital migration, near the mutual 2:1 mean motion resonance (MMR), 
where the period ratio -Psat/-Pjup = 2 (today this ratio is 2.49). The angle Aj^p — 2Asat — 
where w is the perihelion longitude of either Jupiter or Saturn, can then resonate with 
AA, creating widespread chaos around L4 and L5. Small bodies scattered by planets into 
the neighborhood of Jupiter's orbit can chaotically wander near L4 and L5, where they are 
permanently trapped once i-*sat/-fjup moves away from 2. 

A natural consequence of chaotic capture is that orbits fill all available space charac- 
terized by long-term stability. In M05, the planetesimals dynamically evolving from the 
transplanetary disk scatter off of the giant planets, acquire high-inclination orbits, and re- 
main on these orbits after capture. This creates a wide inclination distribution of captured 
bodies, and resolves the long-standing conflict between previous formation theories and ob- 
servations discussed above. 

M05 placed chaotic capture in the context of the original Nice model (hereafter ONM; 
Tsiganis et al. 2005, Gomes et al. 2005; hereafter ONM), where migration of Jupiter 
and Saturn past their mutual 2:1 MMR is thought to trigger an instability during which 
Uranus and Neptune are scattered into the outer planetesimal disk. Their orbits subsequently 
stabilize and circularize near 20 and 30 AU by dynamical friction (Binney & Tremaine 1987). 
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There are many things to hke about the Nice model, but that does not mean it is correct. It 
needs to be continually tested against all possible constraints. Indeed, several inconsistencies 
of the ONM have been already pointed out leading to model's revisions (Morbidelli et al. 

2007, Levison et al. 2011). 

It is now thought that Jupiter and Saturn have not smoothly migrated over the 2:1 MMR. 
Instead, -Psat/-Pjup probably 'jumped' from <2 to >2.3 when Jupiter (and Saturn) scattered 
off of the ice giants (Uranus, Neptune or a similar- mass planet). This model, known as the 
jumping- Jupiter model, is required to explain the secular architecture of the outer planets, 
orbital distribution of asteroids, and dynamical survival of the terrestrial planets (Morbidelli 
et al. 2009b, 2010, Brasser et al. 2009, Minton & Malhotra 2009, Walsh & Morbidelh 2011, 
Agnor & Lin 2012). Moreover, encounters of Jupiter with one of the ice giants are required 
for capture of irregular satellites around Jupiter (Nesvorny et al. 2007). 

M05's chaotic capture does not work in the jumping- Jupiter model, because the reso- 
nances discussed in M05 do not occur, and alternative capture models have not been inves- 
tigated so far. Here we test how Jupiter Trojans can be captured in the jumping- Jupiter 
model. We find that the Trojans are most likely captured immediately after a close encounter 
of Jupiter with an ice giant. As a result of the encounter, ajup changes, sometimes by as 
much as ~0.2 AU in a single jump. This radially displaces Jupiter's L4 and L5, releases the 
existing Trojans, and can lead to capture of new bodies that happen to have semimajor axes 
similar to ajup when the jump occurs. We call this jump capture. 

2. Capture Simulations 

We take advantage of the results published in Nesvorny & Morbidelli (2012; hereafter 
NM12). NM12 performed nearly 10^ numerical integrations of the early Solar System's in- 
stability. The integrations started at the time when the giant planets were already fully 
formed and nebular gas was dispersed (presumably ~3-10 Myr after the birth of the Sun; 
Haisch et al. 2001, Williams & Cieza 2011). At this time, the giant planets were assumed 
to have orbits in mutual MMRs (that have been established during the previous stage of 
convergent planetary migration in the gas disk; Masset & Snellgrove 2001, Pierens & Nelson 

2008, Pierens & Raymond 2011). A disk of planetesimals was placed beyond the outer- 
most ice giant (hereafter transplanetary disk). The dynamical evolution of the planets and 
planetesimals was then tracked, using an A^-body integrator, through and 100-Myr past the 
instability. 

Here we select three cases from NM12. Their properties are illustrated in Figs. [TJ131 In 
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all three cases, the Solar System was assumed to have five planets initially (Jupiter, Saturn 
and three ice giants). This is because NM12 showed that various constraints (such as the 
final orbits of outer planets, survival of the terrestrial planets, etc.) can most easily be 
satisfied when the system starts with five initial planets and one ice giant is ejected during 
the instability (Nesvorny 2011, Batygin et al. 2012). The case with four initial planets 
requires a massive planetesimal disk to avoid losing a planet, but the massive disk also tends 
to produce excessive dynamical damping and long-range residual migration of Jupiter and 
Saturn that violate constraints (Nesvorny 2011, Batygin et al. 2012). It is therefore difficult 
to obtain a plausible planet evolution starting with four planets (NM12 failed to identify any 
in their 2670 four-planet trials). 

A shared property of the selected runs is that Jupiter and Saturn undergo a series 
of planetary encounters with the ejected ice giant. As a result of these encounters, the 
semimajor axes of Jupiter and Saturn evolve in discrete steps. While the semimajor axis 
can an decrease or increase during one encounter, depending on the encounter geometry, the 
general trend is such that Jupiter moves inward, i.e. to shorter periods (by scattering ice 
giant outward), and Saturn moves outward, i.e., to longer periods (by scattering ice giant 
inward). This process leads to just the right kind of Psat/-Pjup evolution during the instability 
(jumping Jupiter; see §1). 

Note that while having the fifth planet is a convenient way to obtain jumping Jupiter, 
the fifth planet is not (in itself) important for capture of Jupiter Trojans. Instead, their 
capture is controlled by the evolution of Jupiter's (and Saturn's) orbit (see §3.1). Therefore, 
if the future studies will identify plausible jumping- Jupiter cases with four planets, it is 
expected that the capture process described here will apply to those cases as well (unless 
the dynamical history of Jupiter's orbit will strongly differ from the one studied here). 

NM12's simulations were performed using the symplcctic integrator known as SyMBA 
(Duncan et al. 1998). The planetesimal disk was resolved by up to 10,000 disk particles in 
NM12, which was sufficient for instabihty calculations, but will be insufficient here where the 
expected capture probabihty is < 10~^ (M05). To deal with this issue, we developed a new 
method that allows us to track the planetary evolution taken from the original SyMBA run, 
and include a very large number of disk particles whose orbits are numerically integrated by 
the swift_rinvs3 code, part of the Swift package (Levison & Duncan 1994). This works as 
follows. 

We first repeat the selected NM12 jobs using SyMBA, and record the planetary orbits at 
1-yr time intervals. Our modified version of swift_rmvs3 then reads the planetary orbits 
from a file, and interpolates them to any required time sub-sampling (generally 0.25 yr, 
which is the integration time step used here in swif t_rmvs3). The interpolation is done 
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in Cartesian coordinates. First, the planets are forward propagated on the ideal Keplerian 
orbits starting from the positions and velocities recorded by SyMBA at the beginning of each 
1-yr interval. Second, the SyMBA position and velocities at the end of each 1-yr interval are 
propagated backward (again on the ideal Keplerian orbits). We then calculate a weighted 
mean of these two Keplerian trajectories for each planet so that progressively more (less) 
weight is given to the backward (forward) trajectory as time approaches the end of the 1-yr 
interval. We verified that this interpolation method produces insignificant errors. 

The swif t_rmvs3 jobs were launched on different CPUs, with each CPU computing the 
orbital evolution of a large number of disk particles (A^^jisk). The initial orbital distribution of 
each particle set was chosen to respect the initial distribution in the original simulation, but 
differed in details (e.g., the initial mean longitudes of particles were random), so that each set 
behaved like an independent statistical sample. This allowed us to build up good statistics. 
To further improve the statistics, particles were cloned upon first reaching the heliocentric 
distance r < 8 AU. Particles reaching r < 8 AU were cloned by adding N^o particles with 
orbits produced by a small (random) perturbation of the orbital velocity vector (relative 
magnitude ~ 10~^). 

We used a modest number of clones in the initial runs (A^cio — 2-4). Upon convincing 
ourselves that a more aggressive cloning leads to correct results (e.g., captured particles come 
from different clones), we used A^isk = 5,000 per CPU, N^o = 19 and a large number of 
CPUs to obtain an effective resolution with 50 million disk particles in Cases 1 and 2, and 
25 million in Case 3 (Table 1). The whole project was concluded over a period of 8 months. 

The numerical integrations described above were run from a few Myr before the in- 
stability to a few Myr after the onset of the instability (total integration timespan of 10 
Myr; Phase 1). A much longer integration was difficult to achieve, because the interpolation 
method described above had large requirements on the computer memory (planetary posi- 
tions saved at 1-yr intervals over 10 Myr represent Gigabytes of data). Moreover, while the 
final planetary orbits obtained in NM12 matched the real orbits pretty well, they differed in 
details. For example, -Psat/-Pjup sometimes ended up being a bit lower (~2.46) than in the 
present Solar System (2.49). This difference, despite being small, would affect the long-term 
stability of Jupiter Trojans (e.g., Robutel & Bodossian 2009). 

We therefore continued the simulations from Phase 1 using a different method. As 
the planets are orbitally decoupled by the end of Phase 1, the integrator does not need to 
deal with the stochastic outcomes of planetary encounters. Instead, the orbital evolution 
of planets during Phase 2 was governed by scattering encounters with disk particles. As 
a result, the planets slowly migrated toward their current semimajor axis values. This 
phase was followed using the swif t_rmvs3 code modified to include forces that mimic radial 
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migration (inward for Jupiter, outward for Saturn). The migration timescale has been set 
to be equal to that in the original SyMBA simulations (~ 30 Myr e- folding timescale). The 
final orbits of Jupiter and Saturn were tuned so that -Psat/-Pjup = 2.49 in the end. Using 
artificial forces we will also slowly damped planetary eccentricities and inclinations (Lee & 
Peale 2002), in a manner that is consistent with the original evolutions. 

Phase 2 simulations were run for 100 Myr. We did not follow all disk particles during 
Phase 2. Instead, we identified Trojan 'candidates' by selecting particles with orbits near 
Jupiter's L4 and L5 at the end of Phase 1. This selection was very liberal in that most 
of the (hundreds of) selected bodies turned out not to be truly stable Trojans. The non- 
selected particles were discarded, which allowed us to cut down the CPU cost of Phase 2. 
One downside was that given the source population of disk particles was removed, no Trojan 
captures could have occurred in our Phase-2 simulations. This should not be a problem, 
however, because the population of particles was already depleted at this stage, and planetary 
evolution during Phase 2 was not favorable for capture (no planetary encounters, no major 
resonance crossings, etc.). 

To test the long-term stability of Trojans surviving at the end of Phase 2, we performed 
an additional numerical integration over 4 Gyr. This Phase-3 integration used the original 
swif t_rmvs3 code and 0.25-yr timestep. We found that the long-term stability requirements 
shaved off about 50% of Trojans from the population that survived at the end of Phase 2. 
The removed particles typically had D > 60°, large e and/or large i. This result is consistent 
with the expected stability of Jupiter Trojans (Levison et al. 1997, Nesvorny & Dones 2002, 
Robutel & Gabern 2006). 

3. Results 

Here we discuss the results of the numerical integrations described in the previous sec- 
tion. The mechanism of jump capture is illustrated in §3.1. We then examine the orbital 
distribution of captured Trojans and compare it with observations (§3.2). The efficiency of 
jump capture and its implications for the size distribution of planetesimals in the transplan- 
etary disk are discussed in §3.3. In §3.4, we point the possible source of asymmetry between 
the populations of L4 and L5 Trojans. 
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3.1. Capture Process 

While the global evolution of planets was similar in the three selected cases (Figs. [1]13]), 
the detailed behavior of Jupiter's orbit differed from case to case (Fig. H]). This difference 
can be important for Trojan capture and is why different cases were selected in the first 
place. In Case 1, several very close encounters between Jupiter and an ice giant occurred 
near the end of the scattering phase (t = 5.736 Myr) producing a cumulative change of ajup 
from 5.53 to 5.2 AU. In Case 2, the scattering phase of Jupiter lasted over 300 kyr with many 
close encounters contributing to changes of ajup. In contrast, Case 3 showed a relatively poor 
history of Jupiter's encounters lasting 40 kyr only. 

By analyzing these different cases we found that capture of most Jupiter Trojans gener- 
ally occurred during the stage of Jupiter's encounters. In Case 1, for example, roughly 70% 
of captured Trojans that survived to the end of Phase 3 (hereafter the stable Trojans) started 
librating around L4 or L5 at t = 5.715-5.736 Myr after the start of Phase 1. This clearly 
coincides in time with the period of Jupiter's encounters with the ice giant (see Case 1 in 
Fig. H]). By analyzing the capture histories in detail we found that ~50% of stable Trojans 
in Case 1 were captured during the closest encounter between Jupiter and the ice giant at 
t = 5.735869 Myr, when ajup jumped from 5.53 to 5.3 AU (the ice giant was scattered to a 
very eccentric orbit result of this encounter). 

The particles captured at t = 5.735869 Myr had special orbits just before the encounter 
(a ~ 5.3 AU, low e and i < 30°). They were scattered to these orbits by previous encounters 
with planets (Fig. |5]). They were subsequently captured in librating trajectories around L4 
or L5 when the Lagrange points got displaced to ~ 5.3 AU as a result Jupiter's semimajor 
axis jump. This is a clear example of jump capture. In addition, roughly 10% of the stable 
Trojans were jump captured during the previous encounter at t = 5.715 Myr when ajup 
increased (Fig. |5]). They survived on librating trajectories during the closest encounter at 
t = 5.735869 Myr only because they had the right libration phase during the encounter (so 
that a ~ 5.3 AU). 

Out of the remaining ~40% of stable Case-1 Trojans roughly 20% were captured after 
the stage of planetary encounters was over {t > 5.8 Myr). Chaotic capture related to 
sweeping over weaker resonances, such as the 7:3, 12:5 and 17:7 MMRs, was responsible for 
these cases. The other ~20% showed a complicated evolution that was difficult to classify. 
These cases probably correspond to jump capture during smaller jumps of ajup, to chaotic 
capture during the irregular evolution of ajup, or to the combination of both. The capture 
statistics in Cases 2 and 3 were broadly similar: 55% of clear jump captures, 5% of clear 
chaotic captures, and 40% unclear in Case 2; and 50% of jump capture, 20% of chaotic 
capture, and 30% unclear in Case 3. 
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3.2. Final Orbits 

Figure [6] shows an example of the orbital history of a stable Trojan after its capture. 
The changes of the libration amplitude seen in panel (a) are related to sweeping resonances 
(Robutel & Bodossian 2009). The libration amplitude stabilized at t > 10 Myr as planetary 
migration slowed down, and the system did not encounter any important resonances when 
-Psat/-Pjup approached 2.49. While the eccentricity can still significantly change after capture 
(panel b), the inclination of captured orbits typically remained nearly constant (panel d). 
This shows that the inclination distribution of Jupiter Trojans is closely related to that of 
planetesimals near 5 AU during the scattering phase. It is wide mainly because of scattering 
encounters of planetesimals with Jupiter and SaturnQ 

The orbital distribution of stable Trojans produced in our simulations very closely 
matches observations (Figures [7] and [8]). The distribution extends down to very small li- 
bration amplitudes, small eccentricities and small inclinations. These orbits are generally 
the most difficult to populate in any capture model. The inclination distribution of captured 
objects is wide, extending up to i ~ 30°, just as needed. 

To carefully compare the inclination distribution obtained in our model with observa- 
tions, we should ideally need to account for the detection efficiency of objects with different 
inclinations. This is because most surveys look near the ecliptic and tend to have lower 
detection efficiencies for orbits with larger inclinations (e.g., Jewitt et al. 2000). We use 
a magnitude cutoff to avoid this problem. According to Szabo et al. (2007), the Trojan 
population should be (nearly) complete up to if ~ 12. Fig. ^ shows an excellent agreement 
between our model inclination distribution and the one obtained with if < 12. 

Interestingly, despite the very different histories of ajup in Cases 1, 2 and 3, discussed in 
§3.1, the orbital distributions of stable Trojans that were obtained in these cases are similar 
(Fig. [7j). This indicates that jump capture is a robust capture mechanism that is expected to 
produce the correct orbital properties of Jupiter Trojans for a wide range of jumping- Jupiter 
evolutions. 

We apphed the Kolmogorov-Smirnov (K-S) test to the orbital distributions in Fig. [8] 
(Press et al. 1992). The K-S test is a statistical measure indicating whether two sets of 
data (the computed and observed orbits of Trojans in our case) are drawn from two different 



^The wide inclination distribution of Jupiter Trojans is therefore unrelated to and cannot be used to con- 
strain the inclination distribution of planetesimals in the original transplanetary disk. In contrast, Neptune 
Trojans suffered smaller inclination perturbations prior to their capture and can be used to this end (e.g., 
Nesvorny & Vokrouhhcky 2009). 
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distribution functions, or whether they are consistent with a single distribution function. 
The K-S test gives 16%, 68% and 63% probabihty that the computed and measured (with 
cutoff H < 12) distributions of D, e and i are statistically the same, respectively. 

The lower K-S probability for D is caused by a modest shortage of model orbits with 
D < 20°. This, however, varies from case to case. In Case 1, the K-S probability for D is 
60%, while it is only 3% in Case 2. The larger probability in Case 1 may be related to the 
fact that most bodies were captured in Case 1 during a single large jump of ajup (§3.1). Such 
a clean jump can more easily produce D < 20° Conversely, in Case 2, a more continuous 
evolution of ajup may not allow enough bodies to evolve to D < 20°. These conclusions will 
need to be checked with better statistics. 



3.3. Capture Efficiency 

We identified 30 stable Trojans in Case 1 (out of N^isk = 5 x 10^ disk particles), 41 in 
Case 2 (A^disk = 5 x 10^), and 17 in Case 3 (A^disk = 2.5 x 10^). This corresponds to the 
mean weighted capture efficiency of P = (7.0 ± 0.7) x 10^'' for each particle in the original 
planetesimal disk (Table 1), where the formal 1 a error was computed assuming the normal 
distributional As there are 25 known Trojans with H < 9 (this sample is complete), this 
indicates that the planetesimal disk contained ~ 25/(7 x 10^^) = 3.6 x 10^ planetesimals 
with H < 9 (corresponding to diameter D = 80 km for a 7% albedo, Grav et al. 2012). 

This is encouraging because it favorably compares with estimates obtained by other 
means (e.g., Nesvorny et al. 2007, Levison et al. 2008, Charnoz et al. 2009, Morbidelli et al. 
2009a). For example, Charnoz et al. (2009) suggested, using the crater record on Saturn's 
moon lapetus (surveyed by the Cassini spacecraft), that the planetesimal disk contained ~10'' 
planetesimals with D > 80 km. Morbidelli et al. (2009), using a synthesis of constraints 
(mainly from the Kuiper belt), proposed that the disk contained ~10^ planetesimals with 
H <9. 

The stable Trojans captured in our simulations sample the full radial extent of the 



two-dimensional K-S test applied to the e-D and i-D distributions obtained in Case 1 gives 25% and 
50% probabilities, respectively, that the computed and measured distributions are the same. 

■^The capture efBciency is given here for each particle in the original transplanetary disk. M05 instead 
reported, quoting, 'capture efhciencies per one particle cycled through the system as the planets migrate 
through unstable Trojan configurations'. M05 found that this corresponds to 3.4 MEarth in the reference 
simulation of Tsiganis et al. (2005). As Tsiganis et al. used Afdisk ^ 35 MEarth, the capture efficiencies 
reported in M05, 1.8 x 10^^ to 2.4 x 10^^, should be divided by '^10 to compare them to our values. 
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transplanetary disk up to ~30 AU (at least we were not able to detect any preferential 
sampling based on our statistics). This shows that the size frequency distribution (SFD) of 
Trojans should be representative for the SFD in the whole transplanetary disk (at least) up 
to ~30 AU (Morbidelh et al. 2009a). 

M05 estimated that the present mass of the Trojan population is Mxro ~ 10~^ Msarth- 
Jewitt et al. (2000), on the other hand, suggested that Mxro ~ 9 x 10^^ MEarth- Using bulk 
density p = 1 g cm~^ (e-g-j Marchis et al. 2006) instead of Jewitt's p = 2 g cm~^, and updated 
albedo (7% from Grav et al. 2012 instead of Jewitt's 4%), we find that MTro ~ 2 xlO"^ MEarth- 
With Mxro ~ (1-2) X 10~^ MEarth, we can therefore estimate that the planetesimal disk mass 
was ~ (1-2) X 10-V(7 X 10"^) = 14-28 MEarth- This is consistent with Mdisk = 20 MEarth 
used in NM12. 



3.4. Asymmetry 

The difference in the number of leading and trailing Trojans is a long-standing problem 
in planetary science (see Marzari & Scholl 2002 for a review). This is because all capture 
mechanisms proposed so far, including chaotic capture of M05, are expected to produce 
symmetric distributions with the numbers of bodies in each swarm, N{L^ and N{L^), being 
nearly the same (up to statistical fluctuations). Planetary migration is also not expected 
to change f^^ = N{Li)/N{L^)^ nor is the long-term stabilityo In addition, the existing 
asymmetry is not simply due to different collisional evolution of the L4 and L5 Trojan 
swarms, because it persists even if known collisional families at L4 are removed (O'Brien & 
Morbidelh 2008). 

Various research groups published estimates of /45 that apply to different limiting mag- 
nitudes/sizes. Szabo et al. (2007) estimated from the Sloan Digital Sky Survey (SDSS) that 
/45 = 1.6 ± 0.1. Survey with the Subaru telescope gave /45 ~ 1.8 (Nakamura & Yoshida 
2008). Estimates from the Wide-field Infrared Survey Explorer (WISE) gave f^^ = 1.4 ±0.2 
(Grav et al. 2011) and f^^ ~ 1.34 for diameters D > 50 km (Grav et al. 2012). Therefore, 
/45 = 1.2-1.8 according to these works. 

Figure [9] shows f^^ of the known Jupiter Trojans as a function of H (data from the Minor 
Planet Center). The ratio is wiggly for H < 10, because only a very few bright Trojans exist. 



''Gomes (1998) claimed that planetary migration can change /45, but these changes were probably related 
to the choice of the initial orbits in Gomes (1998) rather than to an asymmetry of the effects of planetary 
migration itself (e.g., O'Brien 2012). 
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The statistics for if < 9 is apparently not large enough to rule out /is = 1 with confidence. 
At the faint end, on the other hand, the sample is incomplete, and can be influenced by 
a few large coUisions that generated a lot of small debris (e.g., Broz & Rozehnal 2011). 

To assess the significance of the asymmetry we find it best to use the population with 
H < 12 (complete sample according to Szabo et al. 2007). The number of the known 
L4 Trojans with if < 12 is 361, if 9 known Eurybates family members with if < 12 are 
removed (Broz & Rozehnal 2011). For comparison, there are 279 L4 Trojans with H < 12. 
This indicates that = 1.3 ± 0.1, where the formal error was computed as /45[1/A''(L4) + 
l/iV(L5)]°'^. The asymmetry is therefore significant at ~3cr for H < 12. 

The jump capture discussed in §3.1 is potentially capable of producing an asymmetry. 
For example, the asymmetry can arise shortly after planetesimals are jump captured, while 
the ice giant still remains on a Jupiter-crossing orbit. During this time, the ice giant can 
traverse one of the Lagrange swarms and scatter captured bodies, causing the local popula- 
tion to drop. If so, the observed /45 > 1 would indicate that the ice giant traversed the L5 
cloud shortly after the bulk of the Trojan population was captured. Note that this source 
of asymmetry does not apply to chaotic capture described in M05, because the orbits of ice 
giants do not reach down to Jupiter's orbit in M05. 

Prom our simulations, we get f^^ ~ 1.3±0.5 in Case 1, f^A = I//45 = 1.4±0.4 in Case 2, 
and /54 = 1.8 ± 0.9 in Case 3. While all three cases therefore show a formal asymmetry, the 
statistics is not good enough to rule out f^^ = 1 at more than 1 a. It is therefore possible 
that we are just seeing statistical fiuctuations of a small sample. Unfortunately, increasing 
the statistics to 3 cr with the method described in §2 is not computationally feasible at this 
time, because we would need to increase the number of disk particles by a factor of ~10. We 
leave this issue for a future work. 



4. Conclusions 

Here we discussed a new model for capture of Jupiter Trojans. The jump capture, as we 
call it, occurs when Jupiter undergoes a series of scattering encounters with an ice giant, and 
ajup evolves through a number of discrete steps as a result of these encounters. We performed 
numerical integrations of jump capture, where orbits of 125 million disk planetesimals were 
tracked over the period of discrete changes of ajup. The captured bodies were followed from 
the time of their capture, presumably some 4 Gyr ago, to the present time. The number 
and orbits of the surviving bodies were compared with observations of Jupiter Trojans. Our 
results are summarized as follows: 
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1. We found that the efficiency of jump capture is (6-8) x 10~^ for each particle in the 
original transplanetary disk. This, and the number of known Trojans with H <9, im- 
ply that the planetesimal disk contained (3-4) x 10^ bodies with H < 9 (corresponding 
to diameter D = 80 km for a 7% albedo). The inferred mass of the planetesimal disk 

is Mdisk = 14-28 MEarth. 

2. The orbital distribution of stable Trojans obtained in our simulations provides a good 
match to the observed distribution, including orbits with small libration amplitudes, 
small eccentricities and small/large inclinations. The present wide inclination distri- 
bution of Jupiter Trojans reflects the distribution of planetesimals near 5 AU during 
the planetary instability. 

3. The jump capture is potentially capable of explaining the observed asymmetry of 
Jupiter Trojans {N{Li)/N{L^) = 1.3 ± 0.1 for the complete sample with H < 12). 
The asymmetry can be related to (a few significant) late passages of an ice giant near 
L5 that presumably depleted the L5 population. Future modeling work will need to 
improve the capture statistics and test this possibility at a larger statistical confidence 
that it was done here. 

In a broader context, the work presented here provides support for the jumping-Jupiter 
model (Morbidelli et al. 2009b, 2010; Brasser et al. 2009), and shows a good consistency of 
the planetary-instability simulations pubhshed in NM12. 
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Table 1: The statistics of Trojan capture. The rows are: the (1) number of disk particles 
(-^disk), (2) number of captured stable Trojans (iVcap), (3) probability of capture (-Pcap), and 
(4) asymmetry in the number of L4 and L5 Trojans (/45 = N {L4) / N (L^)) . 
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Fig. 1. — Orbital histories of the outer planets in Case 1. The planets were started in the 
(3:2,3:2,2:1,3:2) resonant chain, and Mdisk = 20 MEarth- (a) The semimajor axes (solid lines), 
and perihelion and aphelion distances (dashed lines) of each planet's orbit. The black dashed 
lines show the semimajor axes of planets in the present Solar System, (b) The period ratio 
-Psat/-Pjup- The dashed line shows -Psat/-Pjup = 2.49, corresponding to the period ratio in the 
present Solar System. The shaded area approximately denotes the zone where the secular 
resonances with the terrestrial planets occur. 
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Fig. 2. — Orbital histories of the outer planets in Case 2. See the caption of Fig. [T]for the 
description of orbital parameters shown here. 
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Fig. 3. — Orbital histories of the giant planets in Case 3. See the caption of Fig. [T]for the 
description of orbital parameters shown here. 
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Fig. 4. — Orbital evolution of Jupiter during planetary encounters in Cases 1 (left), 2 (mid- 
dle) and 3 (right). The solid hnes show the semimajor axis (black), and perihelion and 
aphelion distances (gray). The dashed vertical lines delimit the interval of Jupiter's encoun- 
ters with an ice giant. 
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Fig. 5. — Orbital evolution of a disk particle that was captured as a stable Trojan in Case 1: 
(a) semimajor axis (black line), and perihelion and aphelion distances (gray lines) of the 
particle, (b) eccentricity, and (c) inclination. Jupiter's and Saturn's semimajor axes are 
shown in (a) by red and blue lines, respectively. The particle orbit remained near its starting 
location in the transplanetary disk for up to t = 3.5 Myr after the start of the integration. 
The changes of e and i were minor during this stage. Then, at t = 3.5 Myr, the particle 
was scattered by Neptune, evolved inward, and a series of subsequent encounters with ice 
giants raised orbit's e and i to moderate values. At t = 5.4 Myr, particle's eccentricity 
evolved to very high values (e = 0.8) by encounters with Saturn. Finally, shortly before 
t = 5.735869 Myr, when the closest encounter of Jupiter with an ice giant occurred (vertical 
dashed line), the particle was scattered by Jupiter. This changed its orbit in just the right 
way for capture to be possible at t = 5.735869 Myr (i.e., a ~ 5.3 AU and small e prior 
to capture). The large inclination of captured orbit was established by several scattering 
encounters with Jupiter shortly before capture. 
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Fig. 6. — Orbital evolution of the particle shown in Fig. [5] after its capture at t = 
5.712936 Myr (indicated here by the vertical dashed line). The angle AA circulates be- 
fore capture and is not shown in panel (a) for t < 5.712936 Myr. The bright grey line in 
panel (c) is the semimajor axis of Jupiter. 
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Full Libration Amplitude (deg) 

Fig. 7. — Orbits of stable Trojans obtained in Cases 1 (circles), 2 (triangles), and 3 (stars). 
The full libration amplitude D corresponds to the angular distance between extremes of 
^ ~ -^Jup during libration. The black dots show the orbital distribution of real Trojans. 
The proper orbital elements shown here were computed by the method described in Broz & 
Rozehnal (2011). We found no significant difference between the orbital distributions of L4 
and L5 Trojans, and clumped these distributions together. 
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Fig. 8. — The cumulative distribution of the (a) full hbration amplitude, (b) eccentricity, 
and (c) inclination. The various hues shown here denote the model distributions (solid), 
known Trojans (dashed), and known Trojans with absolute magnitude H < 12 (dot-dashed). 
According to Szabo et al. (2007), the population of known Trojans with H < 12 should be 
nearly complete. The difference between the dashed and dot-dashed lines in panel (c) is 
related to the incompleteness of the faint Trojans with high orbital inclinations. 
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Fig. 9. — Asymmetry between the populations of L4 and L5 Trojans. The sohd hnc shows 
/45 = N{Li) /N[L^) for the known Trojans as a function of absolute magnitude H. The 
Eurybates family, as identified by Broz & Rozehnal (2011), was removed from N{L^ (the 
dash-dotted line shows with the Eurybates family included in N[L:i)). The gray region 
denotes our 1 a statistical error estimate of 745. 



